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I review the study of real (Minkowski) time correlators in hot, weakly coupled Yang-Mills theory via lattice 
methods. I concentrate on the Minkowski time topological susceptibility, which is related to the efficiency of 
baryon number violation at high temperature. It can be computed by approximating the IR fields as classical 
and solving their dynamics nonperturbatively on the lattice. However it is essential to include the UV degrees 
of freedom. Their influence can be computed perturbatively and included in the lattice model by the addition of 
auxiliary fields. 



Instanton in QCD 



1. Introduction 

Most of the work to date in lattice gauge theory 
has focused on computing Euclidean quantities. 
While this is a valid pursuit, there are very inter- 
esting questions to be asked about time-like sep- 
arated correlators. For instance, almost all ques- 
tions about dynamics involve such correlators. 
For QCD it appears at this time that we have 
no first principles nonperturbative tools to study 
these questions. However, the situation is some- 
what different in the case of weakly coupled, hot 
Yang-Mills theory. This is a relevant question for 
the SU(2) sector of the standard model, for which 
the coupling is weak, a w ~ 1/30. Naively we 
would think that perturbation theory should then 
work, and nonperturbative tools are not needed. 
However at very high temperatures this is not the 
case. High temperature questions are also inter- 
esting in cosmology; in particular we may need to 
understand the behavior of the SU(2) weak sector 
at high temperatures to understand the origin of 
the cosmological baryon number abundance. 

2. Baryon number violation and topologi- 
cal susceptibility 

First, a big question; how efficiently is baryon 
number violated in the standard model? To see 
that it is violated, recall what happens when there 
is an instanton in QCD: left handed quarks go in, 
right handed quarks go out (see Fig. [l]). 

Instantons are also permitted in SU(2) weak; 
the difference is that only left handed fermions 




Instanton in SU(2) weak 




Figure 1. Top: cartoon of what an instanton 
in QCD does; it turns right handed quarks into 
left handed ones. But because SU(2) weak only 
couples to left handed fermions, a weak instanton 
(Bottom) creates a net particle number. 



couple to SU(2) weak, and so, while left handers 
go out, no right handers go in. This is particle 
number violation; in particular, baryon number 
is violated, as first noted by t'Hooft Q. 

So why don't we all die? Because the decay 
constant for baryons at normal temperatures is 
tiny. In vacuum we can relate the baryon num- 
ber violation rate to the density of instantons in 
the Euclidean theory, giving a decay rate of order 
exp(-167r 2 /.g^) < 10~ 160 . 



2 



The situation may be different at high temper- 
atures. Arnold and McLerran Q have shown that 
the decay lifetime r for a baryon number excess 
is related to the Minkowski space topological sus- 
ceptibility (often called the sphaleron rate), 
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x(^^(^^/3^(0,0)), (1) 

through t- 1 = (39/4T 3 )r. Here T should be 
computed using the thermal density matrix but 
at zero chemical potential. This relation applies 
for temperatures T » m pro ton and baryon num- 
ber densities N B /V < T 3 . 

Note that T is NOT the Euclidean topological 
susceptibility; Eq. (|l|) involves an integral over 
real (Minkowski) time. The two susceptibilities 
may be related in vacuum, but they are certainly 
not at high temperature; the Euclidean suscepti- 
bility is always exponentially small and falls with 
temperature, while T grows with increasing tem- 
perature, eventually becoming 0(a^T 4 ). 

To compute F we need a nonperturbative way 
to measure f(F^ v F^(x,t)F* p F^{0,0)) at large 
T. Because the time integral involves all times, 
we cannot do this reliably by analytic continua- 
tion from Euclidean time. And as we shall now 
see, perturbation theory also proves unreliable. 

3. Thermal Perturbation Theory and Clas- 
sical Physics 

My discussion here will follow ||. One can 
arrange thermal perturbation theory in the real 
time formalism in terms of order symmetrized (F) 
and antisymmetrized (retarded and advanced) 
propagators. Each diagram gives several contri- 
butions in which different propagators are F or 
advanced/retarded, and as we add loops to a dia- 
gram, the maximum possible number of F prop- 
agators increases by 1 per loop. In this formal- 
ism, temperature appears as the replacement on 
F propagators of the vacuum zero point fluctua- 
tion amplitude with the thermal fluctuation am- 
plitude; for bosons, 
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Neglecting interactions E = y/k 2 + m 2 . At high 
temperatures the Higgs condensate, responsible 
for particle masses, dissolves and m ~ 0. In the 
infrared, the Bose distribution function is large 
and we can make the expansion 
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which in the infrared is very well approximated by 
T/E. The same expansion for the fermions starts 
at 0(E/T), so their IR effects are negligible. 

If we start out with the classical theory we get 
the same perturbation theory except that only 
diagrams with the maximum number of F prop- 
agators appear, with Eq. (||) simply replaced by 
T/E. Hence, the IR of a weakly coupled quan- 
tum theory at T > m behaves at leading order 
as a classical theory. 

Moreover, this leads to a failure of perturbation 
theory. In vacuum, the loop counting expansion 
parameter is a Wl which is small; but if one propa- 
gator per loop is an F propagator, the expansion 
parameter becomes (T/E)a w . If some bosonic 
mass scale is m < a w T, then for k ~ a w T, the 
expansion parameter is 1 and perturbation the- 
ory breaks down. However, perturbation theory 
only fails in the IR, and it fails only where the 
fields behave classically up to small corrections. 
Therefore we can treat the strongly coupled IR 
behavior by treating the IR degrees of freedom 
classically. 

4. Classical Theory on the Lattice 

With this in mind, Grigoriev and Rubakov con- 
jectured that any unequal time IR dominated cor- 
relator in hot Yang-Mills theory would take the 
same value as its classical theory analogue, up to 
0(a w ) corrections j|. This turns out to be in- 
correct; it is essential to include the interactions 
between the IR modes and the (quantum but per- 
turbative) UV modes || . But for starters we will 
"throw out" the UV and just look at T in the 
regulated classical theory. 

We need a nonperturbative regulation which 
will preserve exact gauge invariance and remove 
the unwanted (nonclassical) UV degrees of free- 
dom; naturally the lattice is the best candidate. 
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Real time classical Yang-Mills theory on the lat- 
tice looks almost like the molecular dynamics al- 
gorithm applied to the 3-D path integral without 
fermions. To each link matrix ?7, S SU(2) we 
associate a canonical momentum Ei in the Lie 
algebra, satisfying equations of motion || 



dt 
dt 
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which arise from the Hamiltonian 

H = ^Wilson, 3-D + l E i( X ) E i( X ) 



(5) 



(6) 



There are two differences from molecular dyam- 
ics. First, the E fields must satisfy Gauss' Law, 

D L -E = 0, (7) 

which is conserved by the equations of motion 
but must be enforced on the initial conditions and 
considered in the thermalization. Second, we in- 
terpret t not as a "fake" time put in to perform 
a Monte-Carlo, but literally as Minkowski time. 

To determine T we must measure the topology 
of the 4-D configuration gotten by taking 3-D con- 
figurations at successive times as neighboring 3-D 
slices of a 4-D lattice, see Fig. |[ 

/ — — 3-D final configuration 




Time Series 
Built by 
Molecular 
Dynamics 



3-D starting configuration 

Figure 2. The 4-D space for which we must 
measure toplogy. 

The space has open boundary conditions, so 
topology is technically undefined. Also we need 
a very efficient algorithm, since the algorithm 
which generates the configurations is very fast. 



We fix the first problem by repeatedly append- 
ing a "new" initial and final 3-D slice, equal to 
the old one after some 3-D cooling. Sufficient 
cooling leads quickly to the vacuum (because we 
are in 3-D), and topology is defined with vacuum 
boundaries. The topological suceptibility is not 
affected in the limit of large molecular dynam- 
ics time. Finding an efficient algorithm is made 
easier because the gauge configurations are very 
smooth; first, the mean plaquette is much closer 
to 1 than in the QCD setting because of the differ- 
ence between 3-D and 4-D thermodynamics (not 
surprising as the underlying coupling is weak). 
Also, the temporal spacing is very small. Cooling 
within 3-D planes makes topological configura- 
tions spread in the 3 dimensions; they contract in 
the temporal direction, but it is highly resolved 
anyway and this is not a problem. It then re- 
mains to integrate FF with an 0(a 2 ) improved 
operator. The details are in JjJ. 

Is r in the classical theory lattice spacing inde- 
pendent? The data ||, presented in Fig. Q, says 
no. 
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Figure 3. Topological suceptibility of hot lattice 
SU(2) Yang-Mills theory plotted against lattice 
spacing. Linear lattice spacing dependence is ev- 
ident at the smallest lattice spacings. 

This behavior is expected analytically gj; it 
means that the UV must matter, and we should 
treat it properly. 
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5. Putting in the UV physics 

We need to include properly the interactions 
between IR and UV. The UV does not behave 
classically; but it behaves perturbatively, and its 
leading order influence on the IR, the Hard Ther- 
mal Loops (HTL's), is known ||. The HTL ef- 
fective action is nonlocal , but it is equivalent 
to the effect of a large collection of classical, ad- 
joint charged particles pLlj ] . This means we may 
be able to incorporate the effects in a local way, 
and it provides the intuition of why we have to. 

First, the intuition. The UV modes are short 
wavelength and so behave like classical particles. 
A large collection of charged particles is a plasma, 
and plasmas are conducting. By Lenz's Law, con- 
ducting media resist changes in magnetic fields. 
But evolution of IR magnetic fields is precisely 
what is required to get topology change. Hence, 
the more UV modes there are, the smaller F 
should be. This is basis of the arguments of 
Arnold, Son, and Yaffe |Q. 

How should we add the UV effects to the lattice 
IR theory? Two approaches seem reasonable; an 
N-body treatment and a Boltzmann-Vlasov treat- 
ment. 

In the N-body approach, we add a large num- 
ber of "particle" degrees of freedom to the lattice 
gauge theory [|l3|. Each particle has a coordi- 
nate, a momentum, and an adjoint representation 
charge. In Hu and Muller's formulation, the coor- 
dinates and momenta are each continuous degrees 
of freedom, so the particles "live" in the continu- 
ous space between the lattice sites, see Fig. ||. 

The exact rules for the interactions between the 
particles and the lattice degrees of freedom were 
worked out in |l4[ ). Each particle charge resides 
at the nearest lattice point, and the particle inter- 
acts with the lattice fields at the discrete instants 
when it passes from being nearest one site to be- 
ing nearest another. The system has a conserved 
energy and phase space measure and exactly re- 
spects Gauss' Law; and in the small a and small 
particle charge limit it correctly reproduces the 
classical field + HTL effective theory. 

The details of the Boltzmann-Vlasov approach 
appear in Kari Rummukainen's contribution to 
these proceedings Jig], but I briefly review them 
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Figure 4. Cartoon of how "particle" degrees of 
freedom live in the continuous space between lat- 
tice sites. The charge resides at the vertex nearest 
the particle. 

here. Rather than considering particles individ- 
ually, one considers continuous particle popula- 
tion functions f(p,x,q). Linearizing in the de- 
parture from equilibrium, / — > 8f a (p,x), which 
in the current setting can be further simplified 
by integrating over \p\, leaving only its direction 
p = v. The resulting population function is called 
W a (x,v), and up to normalization it represents 
the net charge of particles moving in the v di- 
rection at position x. The Yang-Mills-Boltzmann 
equations in continuum are fl6| 

(D v F^) a = 3a > (8) 

{v,D») ab W\x,v) = v,F^(x), (10) 

where = (l,u) and d£l v ^ is the solid angle 
measure for integrating over v. 

The continuum equations not only need to be 
regulated on a lattice, discretizing space; they 
also need regulation on the sphere, which we do 
by expanding in spherical harmonics and truncat- 
ing the series at a finite ^ max - For even ^ max there 
is surprisingly little truncation error |15|j . 

In either approach, we expect the final answer 
to vary inversely with the conductivity, which is 
proportional to the Debye mass squared m D . It 
makes sense, then, to fit T to the form 

' a 2 T 2 \ 

1 ^<r- (ii) 



■i a 



r = K ' 



with k' to be determined. The factor of (aT) is 
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the four powers of the nonperturbative scale aT 
required by dimensions and the last power of a 
is the effect of the conductivity. The numerical 
data from the two approaches to the UV physics 
are in good agreement, as shown in Figure |5[ 



Figure 5. T, scaled by the expected parametric 
behavior, plotted against reciprocal conductivity. 
Theory predicts a flat line. Data with squares are 
using the N-body approach, triangles are with the 
Boltzmann-Vlasov approach. 



The result is k' ~ 45, which is abundantly fast 
enough to erase any GUT generated initial baryon 
number (unless the GUT violates B — L) and also 
fast enough to make electroweak baryogenesis vi- 
able if there is sufficient CP violation and loss of 
equilibrium at the electroweak phase transition. 

6. conclusion 

Real time techniques are available when the 
coupling is small. Although naively perturbation 
theory should be reliable in this case, it breaks 
down whenever some particle mass comes on or- 
der m ~ aT; but in this case the nonperturbative 
physics occurs only in the infrared, which is clas- 
sical and can be treated nonperturbatively in real 
time on the lattice. 

The influence of UV modes on the IR is very 
important in a gauge theory, and there are now 



two different techniques, an N-body approach and 
a Boltzmann-Vlasov approach, to include it. The 
techniques are in good agreement. 

The Minkowski topological suceptibility of pure 
Yang-Mills theory has been measured accurately, 
and this question is now closed. However there 
are still very interesting questions in nonpertur- 
bative IR physics at the electroweak phase tran- 
sition which remain to be addressed, such as the 
dynamics of (Higgs) scalar condensates during the 
phase transition. 
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